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We present a theory of carrier-induced ferromagnetism in diluted magnetic semiconductors 
(IIIi-^MnxV) which allows for arbitrary itinerant-carrier spin polarization and dynamic correla- 
tions. Both ingredients are essential in identifying the system's elementary excitations and describ- 
ing their properties. We find a branch of collective modes, in addition to the spin waves and Stoner 
continuum which occur in metallic ferromagnets, and predict that the low-temperature spin stiffness 
is independent of the strength of the exchange coupling between magnetic ions and itinerant carriers. 
We discuss the temperature dependence of the magnetization and the heat capacity. 
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Introduction. The recent discovery of carrier-induced 
ferromagnetism in diluted magnetic semiconductors 
(DMS) has generated intense interest, in part because 
it opens the prospect of developing devices which com- 
bine information processing and storage functionalities 
in one material |^|. Critical temperatures T c exceed- 
ing 100K have been realized || by using low-temperature 
molecular beam epitaxy growth to introduce a high con- 
centration c of randomly distributed Mn 2+ ions in GaAs 
systems with a high hole density c*. The tendency to- 
ward ferromagnetism and trends in the observed T c 's 
have been explained using a picture JlO| fTi] in which uni- 
form itinerant-carrier spin polarization mediates a long- 
range ferromagnetic interaction between the Mn 2+ ions 
with spin S = 5/2. 

We present here the first theory which accounts for 
dynamic correlations in the ordered state and is able 
to describe its fundamental properties. We use a path- 
integral formulation and, in the RKKY spirit of the pre- 
vious theory, integrate out the itinerant carriers, and ex- 
pand the effective action for the impurity spins up to 
quadratic order. In addition to the usual Goldstone spin 
waves, we find itinerant-carrier dominated collective exci- 
tations analogous to the optical spin- waves in a ferrimag- 
net. We also find that the spin stiffness is inversely pro- 
portional to the itinerant-carrier mass and independent 
of the strength of the exchange coupling between mag- 
netic ions and itinerant carriers, in contrast to the mean- 
field critical temperature which is proportional to the 
itinerant-carrier mass and to the square of the exchange 
coupling. We extend our calculations to the critical tem- 
perature using an approximate self-consistent spin-wave 
scheme and address the temperature-dependent magne- 
tization and specific heat. 

Hamiltonian and effective action. We study a model 
which provides an accurate description of Mn based 
zincblende DMS's. Magnetic ions with S = 5/2 at po- 
sitions Ri are antiferromagnetically coupled to valence- 
band carriers described by envelope functions, 



H = H a + J pd d 3 rS(f)-s{r) 



(1) 



where S(r) — Si5(r — Rj) is the impurity-spin 
density. The itinerant-carrier spin density is ex- 
pressed in terms of carrier field operators by s(r) — 
\ J2acr' ^K^^CTCT'^CT'Cr)- (t 1S tne vector of Pauli spin 
matrices.) Hq includes the valence-band envelope- 
function Hamiltonian pi and, if an external magnetic 
field B is present, the Zeeman energy, 
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The effective mass, chemical potential, and g-factor of 
the itinerant carriers are labeled by m* , /x*, and g* . The 
model we use here is, thus, related to colossol magnetore- 
sistancc (CMR) materials 17 and identical to those for 
dense Kondo systems, which simplify when the itinerant- 
carrier density c* is much smaller than the magnetic ion 
density c The fact that c* /c <C 1 in ferromagnetic- 
semiconductor materials is essential to their ferromag- 
netism. Similar models have been used for ferromag- 
netism induced by magnetic ions in nearly ferromagnetic 
metals such as palladium [|l9| . 

We represent the impurity spins in terms of Holstein- 
Primakoff (HP) bosons J2(|. By coarse graining, the spin 
density S(r) can be replaced by a smooth function, 



S+(r) = ^2cS-tf(r)b(r)^J b{r) 

S-(r) = tf(r)^2cS -bt(r)b(r) 
S z {r)=cS-b\r)b{r) 



(3) 
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with bosonic fields 6'(r),6(r). The partition function as 
a coherent-state path-integral in imaginary times reads, 



Z= T>[zz]V[W]e 



- j d-rL(zz,*4<) 



(6) 
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with L = Jd 3 r [zd T z + £ CT y a d T V a ] + H(zz, **). The 
bosonic (impurity spins) and fermionic (itinerant carri- 
ers) degrees of freedom are labeled by the complex vari- 
ables z, z and the Grassmann numbers v f', 'J, respectively. 

Since the Hamiltonian is bilinear in fermionic fields, we 
can integrate out the itinerant carriers and arrive at an 
effective description in terms of the localized spin density 
only, Z — J D[zz] exp(—S e g[zz\) with the action 

S e s[zz] = J dr Jd 3 r[zd T z — g/iBB(cS — zz)] 

- lndet [(G*^)- 1 + 5G- X {zz)] . (7) 

Here, we have already split the total kernel G _1 into a 
mean- field part (G MF )~ 1 and a fluctuating part SG~ X , 
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where A = cJ p dS — g* PbB is the zero-temperature spin- 
splitting gap for the itinerant carriers. The physics of the 
itinerant carriers is embedded in the effective action of 
the magnetic ions. It is responsible for the retarded and 
non-local character of the interactions between magnetic 
ions, which is described here for the first time. 

Independent spin-wave theory. The independent spin- 
wave theory, which is a good approximation at low 
temperatures, is obtained by expanding Eq. (Q) up to 
quadratic order in z. We find (in Fourier representation) 



S e s[zz] = — ^2 

\p\<p a ,m 
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in addition to the temperature-dependent mean-field con- 
tribution. A Debyc cutoff (pj? = 6n 2 c) ensures that we 
include the correct number of magnetic ion degrees of 
freedom. The kernel of the quadratic action is the in- 
verse of the spin-wave propagator, 

D^ip, v m ) = -iv m + gp B B + J pd n* 

+ E G f F & "n)Gf IF (k + p, u> n + v m ) (11) 



with the mean-field itinerant carrier Green's function 
G^ IF (k,u n ) = - [iu n - (e^ + aA^- p*)]~\ The 
mean-field spin density is denoted by n* = (n^ — n|)/2, 
and ej: = h 2 k 2 /{2m*). 

Excitations. We obtain the spectral density of the spin- 
fluctuation propagator by analytical continuation, iv m — > 
+ i0+ and A(p,Q) = lmD(p, Q)/tt. We find three 
different types of spin excitations. In all figures we take 
B = and use as typical parameters H m* — 0.5m fi , 



free-electron mass. For these parameters the mean-field 
itinerant-carrier system is fully polarized at T = 0. 

i) Our model has a gapless Goldstone-mode branch 
(see Fig. |l|) reflecting the spontaneous breaking of rota- 
tional symmetry | pl| . Expansion of the T = propa- 
gator at A > ep, where ep is the Fermi energy of the 
majority-spin band, yields for small and large momenta 
the dispersion of the collective modes 
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and = xA(l ~ A/e p ) + C(l/p 4 ). At short wave- 
lengths we obtain the mean-field result xA, the spin- 
splitting of a magnetic ion in the effective field produced 
by fully spin-polarized itinerant carriers. Note that the 
itinerant-carrier and magnetic-ion mean-field spin split- 
tings differ by a factor of x = c*/(2cS) <C 1. At long 
wavelengths the magnon dispersion in an isotropic ferro- 
magnet is proportional the spin stiffness p divided by the 
magnetization M. In the adiabatic limit |l7jj22| ], £j? < A, 
our long- wavelength result reflects a spin-stiffness due en- 
tirely to the increase in kinetic energy of a fully spin- 
polarized band when the orientation has a spatial depen- 
dence, p — c*h 2 /(Am*), and a magnetization M which 
has opposing contributions from magnetic ions and itin- 
erant carriers, M — cS — c*/2 = cS(l — x). In this limit, 
the mean-field critical temperature and the spin stiffness 
have opposite dependences on the itinerant-carrier mass. 

The low-energy mode describes spin waves in the local- 
impurity system. This contrasts with the case of fer- 
romagnetism induced by local moments |l!| where the 
low-energy mode is primarily in the free-carrier part. 
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FIG. 1. Spin-wave dispersion for c* = O.lnm -3 . The short 
wavelength limit is the mean-field result xA. For comparison, 
we show also the result obtained from an RKKY picture. 

ii) We find a continuum of Stoner spin-flip particle-holc 
excitations. They correspond to flipping a single spin in 
the itinerant-carrier system and, therefore, occur at much 
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larger energies ~ A (see Fig. ||). For A > ep and zero 
temperature, all these excitations carry spin S z = +1, 
i.e., increase the spin polarization, and therefore turn 
up at negative frequencies in the boson propagator we 
study. (When A < ep, excitations with both 5 2 = +1 
and 5 2 = — 1 contribute to the spectral function.) This 
continuum lies between the curves —A — e p ±2 y j'e p ep and 
for A < also between — A + e p ± 2y 1 e p (ep — A). 
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FIG. 2. Stoner excitations and collective modes in the 
free-carrier system for c* = 0.1nm~ 3 . 

iii) We find additional collective modes associated pri- 
marily with the itinerant-carrier system at energies below 
the Stoner continuum (see Fig. 0). At T = we obtain 
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for A > ep. The spectral weight of these modes is 
—x/(l — x) at zero momentum. 

The finite spectral weight at negative energies indicates 
that the ground state is not fully spin polarized because 
of quantum fluctuations. 

Comparison to RKKY picture. For comparison we 
evaluate the T = magnon dispersion assuming an 
RKKY interaction between magnetic ions. This approxi- 
mation results from our theory if we neglect the spin po- 
larization in the itinerant carriers and evaluate the static 
limit of the resulting spin-wave propagator. The Stoner 
excitations and magnons shown in Fig. |^ do not emerge 
and the spin-wave dispersion (Fig. ^) is incorrect except 
for the limit A <C ep. While the RKKY picture does, in 
some circumstances, provide a realistic estimate of T c , it 
completely fails as a theory of the ferromagnetic state. 

Comparison to a ferrimagnet. Some features of the 
excitation spectrum are, not coincidentally, like those of 
a localized spin ferrimagnet with antiferromagnetically 
coupled large spin 5 (corresponding to the magnetic ions) 
and small spin s (corresponding to the free carriers) sub- 
systems on a bipartite cubic lattice. Representing the 



spins by HP bosons and expanding up to quadratic order, 
we either diagonalize the resulting Hamiltonian directly 
or, as in our ferromagnetic-semiconductor calculations, 
integrate out the smaller spins using a path-integral for- 
mulation to obtain equivalent results. We find that 
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with x = s/S and A = 6J5, where J is the exchange 
coupling, and jp~ — (1/3) ^Jl — cos(pia)] with lattice 
constant a. We recover two collective modes, the coupled 
spin waves of the two subsystems. One is gapless, the 
other one gapped with A(l — x). The bandwidth is ccA, 
and the spectral weights at zero momentum are 1/(1 — x) 
and —x/(l — x), respectively, as in our model. 

Self- consistent scheme. Near the transition temper- 
ature T c , the spin-wave density is of the order of c5, 
and expanding around a fully-polarized state is not a 
good starting point. Instead, we adopt a scheme which 
accounts self-consistently for the reduction of spin den- 
sity in the localized impurity and itinerant carrier sub- 
systems [^3|. Our approach is motivated by the Weiss 
mean-field theory, which is expected to be accurate for 
a model with static long-range interactions between the 
magnetic ions. This model can be obtained in our ap- 
proach by taking the Ising limit (i.e., replacing 5 ■ s by 
S z s z ) and letting the mean-field itinerant-carrier spin 
splitting, A(T) = J pc t(S z ), reflect the thermal suppres- 
sion of the impurity-spin density (S z ). In this simple 
model the spin-wave spectrum is independent of momen- 
tum, Q P (T) — J p dn*(T), allowing a unitary transforma- 
tion to spin-wave eigenstates with a site label. The con- 
straint on the number of spin bosons (< 25) on a site is 
then easily applied and we obtain 
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Ip1<Pc 

where Bs(x) is the Brillouin function, and n(u>) is the 
Bose function. The first two terms in the second form of 
Eq. (|l^) give the spin density obtained by treating spin 
waves as bosons, while the last term is the correction 
from spin kinematics which rules out unphysical states. 

Our self-consistent spin-wave approximation consists 
of using Eq. ( |l5| ) when fl p is p— dependent (see Fig. ||). 
Although the characteristic T 3//2 law for the localized- 
ions magnetization is recovered at low temperatures, the 
prefactor is reduced by 1 — (25+ 1) -1 / 2 , compared to the 
correct [ p4| value of linearized spin- wave theory. That is, 
constraining boson populations in momentum space is 
too restrictive at low temperatures, although the error 
is not serious for 5 = 5/2. The transition temperature 
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of the self-consistent spin-wave theory is lower than the 
estimate from Weiss mean-field theory. 

The presence of spin waves also shows up in the specific 
heat CV, which we calculate within the same scheme (see 
Fig. |^) . Calculating the entropy S we obtain 
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where N p = S — SBs(f3S£l p ) is the average number 
of spin bosons at each momentum. The specific heat 
from the kinetic energy of itinerant carriers turns out to 
be negligibly small. The specific heat of magnetic ions 
Cy is proportional to T 3 / 2 at low temperature, shows a 
Schottky-like anomaly, and has a jump at the transition 
temperature. We estimate that the jump is 5% of the lat- 
tice specific heat when T c ~ 100K, suggesting that this 
should be observable. 




FIG. 3. Average impurity and itinerant-carrier spin as a 
function of temperature for c* = 0.1nm~ 3 in Weiss mean-field 
(Weiss) and self-consistent spin-wave theory (SSW). Inset: T c 
as a function of c*. 
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